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ABSTRACT 

High-latitude landscape evolution processes 
have the potential to preserve old, relict sur- 
faces through burial by cold-based, nonero- 
sive glacial ice. To investigate landscape his- 
tory and age in the high Arctic, we analyzed 
in situ cosmogenic 1(l Be and 26 A1 in 33 rocks 
from Upernavik, northwest Greenland. We 
sampled adjacent bedrock-boulder pairs 
along a 100 km transect at elevations up to 
1000 m above sea level. Bedrock samples gave 
significantly older apparent exposure ages 
than corresponding boulder samples, and 
minimum limiting ages increased with eleva- 
tion. Two-isotope calculations ( 26 Al/ ,0 Be) on 20 
of the 33 samples yielded minimum limiting 
exposure durations up to 112 k.y., minimum 
limiting burial durations up to 900 k.y., and 
minimum limiting total histories up to 990 k.y. 
The prevalence of “Be and 26 A1 inherited from 
previous periods of exposure, especially in 
bedrock samples at high elevation, indicates 
that these areas record long and complex sur- 
face exposure histories, including significant 
periods of burial with little subglacial erosion. 
The long total histories suggest that these high- 
elevation surfaces were largely preserved be- 
neath cold-based, nonerosive ice or snowfields 
for at least the latter half of the Quaternary. 
Because of high concentrations of inherited 
nuclides, only the six youngest boulder sam- 
ples appear to record the timing of ice retreat. 
These six samples suggest deglaciation of the 
Upernavik coast at 11.3 ± 0.5 ka (average ± 
1 standard deviation). There is no difference 
in deglaciation age along the 100 km sample 
transect, indicating that the ice-marginal posi- 
tion retreated rapidly at rates of -120 m yr _1 . 


INTRODUCTION 

The surface morphology of many high- 
latitude landscapes is a product of not one, but 
numerous glacial-interglacial cycles (Sugden, 
1977, 1978; Sugden and Watts, 1977). Because 
areas in Greenland, the Canadian Arctic, and 
Scandinavia are periodically covered by cold- 
based glacial ice and snowfields, bedrock sur- 
faces evolve over hundreds of thousands of 
years and through repeated burial by nonerosive 
glacial ice (Anderson et al., 2000; Bierman et 
al„ 1999; Blake, 1970; Briner et al., 2003, 2006; 
Harbor et al., 2006; Kleman, 1992; Kleman and 
Borgstrom, 1994; Marquette et al., 2004; Stro- 
even et al., 2002a, 2002b; Sugden et al., 2005). 
This implies that high-latitude landscapes may 
preserve relict surfaces that formed hundreds of 
thousands or even millions of years ago. 

Study of these old glacial landscapes is 
important for understanding the time scales and 
processes of high-latitude landscape develop- 
ment; however, such studies are challenging 
because the complex exposure and burial his- 
tories of rock surfaces in these areas violate the 
assumptions of many chronological techniques 
such as traditional 10 Be surface exposure dating. 
Here, we employ analysis of paired cosmogenic 
nuclides (“Be and 26 A1) in adjacent bedrock 
and boulder samples to study the age and his- 
tory of the landscape in Upemavik, northwest 
Greenland. We aim to understand the age of 
rock surfaces, durations of exposure and burial, 
efficiency of subglacial erosion, and rate of 
latest Pleistocene ice-margin retreat in a land- 
scape characterized by deep fjords and exten- 
sive uplands. 

Glaciers can erode their beds in a number of 
different ways, including abrasion and pluck- 


ing. However, most forms of subglacial ero- 
sion require basal ice to be at its melting point; 
this enables processes such as basal sliding, 
regelation, and freeze-on of rock and sediment 
to operate (Herman et al., 2011). Warm-based 
glaciers, for which the basal temperature is at 
the pressure melting point, can erode their beds, 
but under cold-based glaciers, where the basal 
temperature does not reach the pressure melt- 
ing point, little erosion occurs (Sugden, 1978). 
Although exceptions have been noted (Atkins 
et al., 2002; Cuffey et al., 2000; Waller et al., 
2012), it is generally accepted that episodes of 
burial caused by cold-based glacial ice are inca- 
pable of removing sufficient material to expose 
fresh rock surfaces after glaciation. 

Measurement of in situ-produced cosmo- 
genic nuclides can provide information about 
landscape history and chronology (Balco, 
2011; Fabel and Harbor, 1999). Production of 
“Be and 26 A1 is dominated by nuclear spall- 
ation, which occurs when high-energy par- 
ticles impact the nuclei of 16 0 and 28 Si atoms 
in quartz, respectively. This impact causes 
material to be ejected from the nucleus, leav- 
ing behind the lighter isotopes, “Be and 26 Al. 
Spallation-induced cosmogenic nuclide pro- 
duction in rock surfaces decreases exponen- 
tially with depth, with the rate of decrease 
depending on the density of the material. The 
attenuation length of high-energy neutrons is 
-160 g cm -2 ; thus, the highest concentrations 
of cosmogenic nuclides are produced in the 
upper several meters of rock. Efficient glacial 
erosion removes this preexisting inventory by 
abrading or plucking the rock’s surface. 

In areas where glacial ice has not completely 
eroded the preexisting exposure history, analy- 
sis of paired “Be and 26 Al data can be used to 
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make inferences about landscape age, exposure 
and burial durations, and subglacial erosion 
efficiency (Bierman et al., 1999; Briner et al., 
2003;, 2006; Briner and Swanson, 1998; Goeh- 
ring et al., 2010; Harbor et al., 2006; Marquette 
et al., 2004; Nishiizumi et al., 1991; Roberts et 
al., 2009; Stone and Ballantyne, 2006; Stroeven 
et al., 2002b; Sugden et al., 2005). If a land- 
scape also contains areas where glacial ice has 
stripped at least -2-3 m of material from rock 
surfaces (Balco, 2011), sufficient to remove 
most nuclides from previous periods of expo- 
sure, 10 Be surface exposure dating of glacially 
transported boulders and abraded bedrock can 
give estimates of the timing and rate of ice- 
margin retreat (Briner et al., 2009; Gosse et al., 
1995a; Hughes et al., 2012; Kelly et al., 2008; 
Nishiizumi et al., 1991; Phillips et al., 1990; 
Young et al., 2011). 

To understand the history of the landscape 
near Upernavik, northwest Greenland, we used 
an experimental design incorporating paired 
bedrock and boulder samples, paired 10 Be and 
26 Al analyses, and samples collected at differ- 
ent elevations and different distances from the 
ice-sheet margin. This approach allowed us to 
make inferences about both subglacial erosion 
efficiency and the timing of ice-margin retreat 
during the latest Pleistocene deglaciation. 

STUDY SITE AND PREVIOUS WORK 

The Upernavik region (~73°N, 55°W; Fig. 1) 
is characterized by dissected terrain that encom- 
passes table-top-shaped islands up to 1 km 
above sea level separated by deep fjords, includ- 
ing the large Upemaviks Isfjord to the north. 
Fjord depths are not well constrained; however, 
water depth in Upemaviks Isfjord is thought to 
be -700 m (Peterson, 2003). The present-day 
ice margin lies -100 km to the southeast of the 
outermost islands. The bedrock in the Uper- 
navik area is predominately granite and granite 
pegmatite, some of which is weakly foliated; 
younger basalt flows cover the granite toward 
the southeastern end of the field area (Escher 
and Pulvertaft, 1995). 

A pronounced contrast exists between the 
low-elevation and high-elevation landscapes. 
Much of the low-elevation bedrock is glacially 
streamlined (Fig. 2), although striations and 
glacial polish were not observed at the loca- 
tions we visited. Conversely, many of the high- 
elevation bedrock forms exhibit pronounced 
surface weathering features (Fig. 3), including 
pedestals up to 10 cm in height, weathering pits, 
and pervasive frost shattering. At high eleva- 
tions, much of the bedrock erodes by breaking 
into sheets that range from several centimeters 
to several tens of centimeters thick, and in many 


cases surfaces are mantled with grus. Both low- 
and high-elevation surfaces are covered with 
perched boulders. 

Previous dating constraints on deglaciation of 
the Upernavik region are based on radiocarbon 
measurements in detrital organic material (Fig. 
1); however, data from this region are relatively 
sparse in relation to the rest of Greenland’s coast 
(Bennike and Bjorck, 2002). A basal radio- 
carbon age from a lake on Svartenhuk Halvp, 
the large peninsula of land -1° south of Uper- 
navik, provides a minimum deglaciation age of 
10.4-10.2 cal ka B.P. (Bennike, 2000). Farther 
away in Melville Bugt, -2° north of Upernavik, 
a whale vertebra found on recently deglaciated 
terrain yielded a radiocarbon age of ca. 9.3- 
9.0 cal ka B.P. (Bennike, 2008). The vertebra 
was likely deposited on the ocean floor after 
the whale’s death and incorporated into the ice 
sheet during a later readvance. Bennike (2008) 
infers that, by ca. 9. 3-9.0 ka, the ice-sheet mar- 
gin at Melville Bugt was similar to or behind 


its present-day position, indicating that the 
coast had been deglaciated by the early Holo- 
cene. Also near Melville Bugt, a basal radio- 
carbon age from a lake sediment core yielded 
a minimum deglaciation age of ca. 9.9-9. 3 cal 
ka B.P. (Bennike, 2008; Fredskild, 1985), and a 
radiocarbon age from shell fragments in marine 
sediment gave an age of ca. 9. 0-8. 5 cal ka B.P. 
(Bennike, 2008; Kelly, 1980). 

USING COSMOGENIC NUCLIDES TO 
UNDERSTAND GLACIAL EXPOSURE 
AND EROSION 

Theoretical Background 

Measurement of in situ-produced cosmo- 
genic nuclides is a useful technique for under- 
standing the timing of surface exposure and 
burial, as well as the efficiency of subglacial 
erosion. Because cosmogenic radionuclides, 
such as 10 Be and 26 Al, form at known rates in 



Figure 1. Regional map of central-western and northwestern 
Greenland. The small box shows the location of the study site as 
shown in Figure 6. Black dots indicate radiocarbon ages as de- 
scribed in the text. 
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Figure 2. Low-elevation locations are characterized by streamlined bedrock forms and 
abundant perched boulders. No glacial striations or polish were observed. Photograph by 
P. Bierman. 


minerals bombarded by cosmic rays (Lai, 
1988), quantifying the nuclide concentration 
in a sample allows inferences to be made about 
that sample’s exposure history (Balco, 2011; 
Gosse and Phillips, 2001). In the simplest case, 
a surface is rapidly exposed to cosmic radiation 


after being buried beneath material sufficiently 
thick to block most cosmogenic nuclide pro- 
duction (e.g., many meters of glacial ice). The 
concentration of cosmogenic nuclides can then 
be used to determine when exposure occurred if 
the surface has not been eroded (Balco, 2011). 


Analysis of two cosmogenic isotopes with dif- 
ferent half-lives can provide information about 
complex exposure and burial histories (Granger 
and Muzikar, 2001; Nishiizumi et al., 1989, 
1991 ). In the case of the two isotopes employed 
here, 26 Al is produced -6.75 times as rapidly as 
10 Be, although the uncertainty of the production 
ratio is not well constrained, and recent work 
has suggested that the ratio may actually be 
closer to -6.9-7.0 (Argento et al., 2013). Once a 
sample is buried and nuclide production ceases, 
the 26 Al/ 10 Be ratio will drop from 6.75 because 
the half-life of 26 Al (0.71 m.y.) is less than that 
of 10 Be (1.36 m.y.). If a sample is exposed again 
following burial, production resumes, and the 
26 Al/ 10 Be ratio increases. 

The 26 Al/ 10 Be ratio data are commonly shown 
on a two-isotope plot, in which the vertical axis 
is the isotope ratio and the horizontal axis is the 
concentration of one of the measured isotopes, 
usually 10 Be (Fig. 4). On this diagram, samples 
that have been continuously exposed plot along 
a “constant exposure” line characterized by 
slowly decreasing 26 Al/ 10 Be ratio but increas- 
ing 10 Be inventory (Fig. 4). When a sample is 
buried, its position on the diagram travels left- 
ward and downward from the constant exposure 
line, along a path characterized by both decreas- 
ing 26 Al/ 10 Be ratio and 10 Be inventory (Fig. 4). 
Although any inverse solution is nonunique (the 
measured isotope concentrations could result 
from an infinite number of different exposure- 
burial scenarios), two-isotope data can be used 
to calculate the most parsimonious solution: a 
minimum duration of initial exposure followed 
by a minimum duration of burial (Bierman et 
al., 1999; Briner et al., 2006; Hakansson et al., 
2008; Lilly et al., 2010; Stroeven et al., 2002b; 
Sugden et al., 2005). 

If samples have been re-exposed follow- 
ing burial (for example, during the Holocene 
period), the assumption of a single period of 
exposure followed by a single period of burial 
is violated. In these cases, minimum limiting 
exposure durations, burial durations, and total 
histories can be corrected for Holocene expo- 
sure by excluding the most recent period of 
nuclide production from the calculations. The 
minimum age of perched boulders free of inher- 
ited nuclides can provide an estimate of the tim- 
ing of deglaciation and can be used to determine 
the duration of Holocene exposure for which 
this correction should be applied. 

Application of the Two-Isotope Approach 
to Glacial Landscapes 

The two-isotope approach has proven to be 
useful for investigating landscapes where sur- 
faces at high elevation appear to be significantly 


iff* 



Figure 3. High-elevation locations exhibit evident surface weathering features including 
pedestals, weathering pits, frost shattering, abundant grus, and sheet weathering. Photo- 
graph by P. Bierman. 
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Figure 4. Schematic depiction of how two-isotope cosmogenic data are plotted, with explanations about important 
processes and the ways in which they are manifested in the isotope data. 


more weathered than surfaces at low eleva- 
tion (Dyke, 1979; Sugden, 1977, 1978). In 
particular, this technique can be employed to 
distinguish between two different scenarios: 
continuous exposure of uplands to subaerial 
weathering versus long-term preservation of 
uplands beneath intermittent cover by nonero- 
sive glacial ice. 

If the ice is not thick enough to cover the land- 
scape completely, the highest upland surfaces 
are continuously exposed as nunataks (Dyke, 
1979). In this case, high-elevation samples have 
old single-nuclide minimum limiting exposure 
ages, close agreement between bedrock and 
boulder samples, and 26 Al/ 10 Be ratios that are 
indistinguishable from constant exposure. This 
scenario has been documented in southwestern 
Greenland (Roberts et al., 2009); such high- 
elevation surfaces protruded through the ice 
sheet, allowing inferences to be made about ice- 
surface elevation during the last glacial period 
(Hakansson et al., 2007; Rinterknecht et al., 
2009; Roberts et al., 2008, 2009). 

Conversely, elevation-dependent weathering 
differences can also result from burial of high- 
elevation surfaces by nonerosive, cold-based ice 
(Fig. 5), wherein highlands are covered by gla- 
cial ice for extended periods of time, but the ice 
is frozen to the bed and cannot perform erosion, 
leaving the landscape unaltered (Sugden, 1977, 
1978; Sugden and Watts, 1977). This scenario 
could occur when ice cover over the lowlands is 
thick (allowing the base of the ice sheet to reach 
its pressure melting point), but ice cover over 
the highlands is thin (prohibiting the base of 
the ice sheet from reaching its pressure melting 


point). In this case, high-elevation samples have 
old single-nuclide minimum limiting exposure 
ages, poor agreement between bedrock and 
boulder samples, and 26 Al/ 10 Be ratios signifi- 
cantly below the constant exposure ratio that are 
indicative of burial (Fig. 5). Isotopic evidence 
for elevation-based erosion discrepancies has 
been found in Baffin Island, Canada (Bierman et 
al., 1999; Briner et al., 2003, 2006), Newfound- 
land, Canada (Gosse et al., 1993, 1995b), and 
Scotland (Phillips et al., 2006). 

METHODS 

Experimental Design 

In the Upemavik region, we collected 20 
bedrock and 1 3 boulder samples for analysis of 
in situ cosmogenic 10 Be and 26 Al. Samples were 
collected along a 100 km northwest-to-south- 
east transect extending from a nunatak -2 km 
inside the present-day ice margin to the farthest 
outboard island (Table 1; Fig. 6). This transect 
is parallel to large fjords immediately to the 
north (Upernaviks Isfjord) and south (Lakse- 
fjorden) and roughly approximates paleo-ice- 
flow direction. Our sampling scheme, known 
as “dipstick sampling” (Stone et al., 2003), 
involved collecting samples at a variety of 
elevations at several locations along the tran- 
sect normal to the ice margin. In Upernavik, 
we collected paired bedrock and boulder sam- 
ples from a range of elevations (-20-1000 m 
above sea level [a.s.l.]) at each “dipstick.” 
Such a sample distribution yields information 
about both vertical and horizontal variability in 


nuclide concentrations as well as differences 
between the exposure history of bedrock and 
overlying boulders, allowing us to make infer- 
ences about ice extent, subglacial erosion, and 
exposure history. 

We present five different types of ages calcu- 
lated with cosmogenic nuclides, all of which are 
minimum limits. For single-nuclide interpreta- 
tions, we present both 10 Be and 26 Al ages, which 
are concordant in simple exposure scenarios. 
For paired-nuclide interpretations, we present 
modeled exposure durations, burial durations, 
and total history (the sum of exposure and 
burial durations). We correct the paired nuclide 
interpretations for Holocene exposure after the 
last deglaciation. 

Field Methods 

Wherever possible, we sampled bedrock 
and boulder pairs in close proximity, usually 
less than 10 m apart. We collected the upper- 
most several centimeters of material from the 
top surfaces of bedrock forms and perched 
boulders using a chisel and hammer (Table 1). 
We recorded latitude and longitude and eleva- 
tion data with a handheld Garmin 12 global 
positioning system (GPS) that has a positional 
uncertainty of <10 m and an elevation uncer- 
tainty of <25 m. 

Laboratory Methods 

Samples were prepared as outlined in Corbett 
et al. (2011). We isolated quartz from the rock 
using a series of both physical and chemical 
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Figure 5. Illustration of bedrock-boulder minimum limiting age relationships and possible scenarios to explain them. 


processes, and Be and A1 from the quartz through 
a series of chemical processes. We used two dif- 
ferent 9 Be carriers in the preparation of these sam- 
ples: a commercial (SPEX) carrier for the first six 
samples and a low “Be/’Be beryl carrier made at 


the University of Vermont for all subsequent sam- 
ples. We added -250 pg of 9 Be to each sample as 
a spike (see GSA Data Repository 1 ). 

In order to ensure that enough A1 was pres- 
ent in the sample for isotopic analysis (2000 pg). 


we added commercial A1 carrier (SPEX) when 
necessary. Immediately after quartz dissolu- 
tion, we removed two aliquots representing 
-6% of the material from each sample. We 
used these replicate aliquots to quantify the A1 


'GSA Data Repository item 2013264, additional methodological information about sample isotopic analysis and retreat rate simulations, is available at http://www 
.geosociety.org/pubs/ft2013.htm or by request to editing@geosociety.org. 
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TABLE 1. SAMPLE LOCATION INFORMATION AND FIELD OBSERVATIONS 


Dipstick number 

Sample name 

Sample type 

Relative elevation* 

Elevation 
(m a.s.l. )t 

Latitude 

(°N)t 

Longitude 

(°W)t 

Rock type 

Thickness 

(cm) 5 

Dipstick 1 

GU091 

Bedrock 

Low 

24 

72.78189 

56.58601 

Granite 

2 

Offshore island 

GU092 

Boulder 

Low 

21 

72.78199 

56.58644 

Granite 

1.5 

Dipstick 2 
Offshore island 

GU093 

Bedrock 

High 

91 

72.73955 

56.38520 

Granite 

2.5 

GU094 

Boulder 

High 

90 

72.74016 

56.38639 

Granite 

2 

GU095 

Bedrock 

Low 

20 

72.74378 

56.38154 

Granite 

3 


GU096 

Bedrock 

High 

778 

72.75386 

55.87299 

Granite 

3 


GU097 

Boulder 

High 

774 

72.75388 

55.87236 

Granite 

2 

Dipstick 3 

GU098 

Bedrock 

Mid 

372 

72.79260 

55.93168 

Granite 

3 

Fjord-dissected terrain 

GU099 

Boulder 

Mid 

366 

72.79308 

55.93311 

Granite 

1 

GU100 

Bedrock 

Low 

74 

72.81168 

55.82623 

Granite 

1 


GU101 

Boulder 

Low 

71 

72.81164 

55.82584 

Granite 

3 


GU102 

Boulder 

High 

980 

72.71763 

55.47279 

Foliated granite 

2 


GU103 

Bedrock 

High 

998 

72.71811 

55.47547 

Granite 

2 


GU106 

Bedrock 

High-mid 

498 

72.74086 

55.55160 

Granite 

4 

Dipstick 4 

GU107 

Boulder 

High-mid 

500 

72.74085 

55.55119 

Foliated granite 

2 

Fjord-dissected terrain 

GU104 

Boulder 

Low-mid 

253 

72.77333 

55.50288 

Granite 

2 

GU105 

Bedrock 

Low-mid 

270 

72.77360 

55.50378 

Granite 

1 


GU108 

Bedrock 

Low 

33 

72.78111 

55.44319 

Granite 

2 


GU109 

Boulder 

Low 

27 

72.78024 

55.44223 

Granite 

2 


GU1 10 

Bedrock 

High 

745 

72.66094 

55.12159 

Granite 

1 

Dipstick 5 

Fjord-dissected terrain 

GU1 11 

Bedrock 

Mid 

325 

72.68163 

55.02502 

Granite 

1 

GU1 12 

Boulder 

Mid 

325 

72.68163 

55.02502 

Granite 

1.5 

GU1 13 

Bedrock 

Low 

90 

72.65953 

54.98345 

Granite 

1 


GU1 14 

Boulder 

Low 

91 

72.65953 

54.98345 

Granodiorite 

1.5 


GU001 

Bedrock 

High 

603 

72.53587 

53.73338 

Granite 

3 

Dipstick 6 

GU002 

Boulder 

High 

603 

72.53587 

53.73338 

Granite 

3 

Ice-marginal terrain 

GU006 

Bedrock 

Low 

539 

72.53920 

53.73214 

Granite 

3.5 

GU017 

Boulder 

Low 

539 

72.53926 

53.73214 

Granite 

1 


GU041 

Bedrock 

High 

898 

72.61509 

53.58872 

Granite 

4 

Dipstick 7 

GU042 

Bedrock 

High-mid 

895 

72.61546 

53.59240 

Granite 

7 

GU043 

Bedrock 

Mid 

857 

72.61594 

53.59608 

Granite 

3 


GU044 

Bedrock 

Low-mid 

808 

72.61743 

53.59646 

Granite 

3 


GU045 

Bedrock 

Low 

776 

72.61870 

53.59650 

Granite 

3 

•Refers to the relative elevation of samples within each vertically oriented “dipstick.” 

f Elevations and locations were recorded in the field with a Garmin global positioning system (GPS) 12, in m above sea level (a.s.l.). 
§ Refers to the thickness of the slab that was collected for analysis. 
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Figure 6. Study site map. Bedrock and boulder samples were collected along a 100-km-long northwest-southeast 
transect near Upernavik, Greenland, at seven different “dipsticks.” Wherever possible, bedrock-boulder pairs 
were taken from a variety of elevations shown by gray ovals. Within each sample pair, the bedrock sample is listed 
first, followed by the boulder sample. In the few instances where boulders were not available, only the bedrock 
sample is listed. Satellite image is from the U.S. Geological Survey. 
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concentration of each sample through induc- 
tively coupled plasma-optical emission spec- 
trometry (ICP-OES) using internal standard- 
ization (with Ga and Y) and measuring two 
different emission lines. The two measurements 
of sample A1 concentration derived from these 
two aliquots were usually within 1%, and always 
within 3.5%, of each other. Measured A1 concen- 
trations of blanks were within 1% of expected A1 
concentrations based on the mass of 27 A1 carrier 
that was added. Thus, the uncertainty in the 26 A1 
concentrations we report predominately reflects 
statistical limitations of accelerator mass spec- 
trometry (AMS) measurements. 

Isotopic Analysis 

We measured isotopic “Be/’Be and “Al/^Al 
ratios by AMS at Lawrence Livermore National 
Laboratory (Perren et al., 2012; Rood et al., 2010). 
We analyzed all 33 samples for 10 Be, and normal- 
ized the samples to standard 07KNSTD3110, 
with a reported “Be^Be ratio of 2850 x 10~ 15 
(Nishiizumi et al., 2007). Measured ‘“Be^Be sam- 
ple ratios ranged from 2.6 x 10 -14 to 1.0 x 1(L 12 
(see GSA Data Repository 1 ), and 10 Be AMS mea- 
surement precisions, including propagated blank 
uncertainties, ranged from 1.9% to 5.4% (let). We 
subtracted process blank ratios (see GSA Data 
Repository 1 ) from sample ratios and propagated 
blank uncertainties in quadrature. 

After determining 10 Be concentrations and 
modeling minimum limiting exposure ages, we 
chose 20 samples to analyze for 26 Al. Sixteen 
of these 20 samples were chosen because of 
their old 10 Be ages (>20 ka); for these samples, 
26 Al analysis aids in distinguishing continuous 
exposure from complex exposure. Independent 
chronology from radiocarbon dating suggests 
that the Upernavik region was deglaciated ca. 
10 cal k.a. B.P. (see “Study Site and Previous 
Work” section); therefore, we inferred that sam- 
ples with 10 Be ages appreciably above ca. 10 ka 
likely do not have simple exposure histories. 
We also analyzed 26 Al in the two youngest bed- 
rock samples (GU001 and GU098) and the two 
youngest boulder samples (GU002 and GUI 04); 
these four samples serve as a validation of sim- 
ple exposure history in young samples and are 
herein referred to as “control” samples. 

For 26 Al analysis, we normalized all samples 
to standard KNSTD10650, with a reported 
26 A1/ 27 A1 ratio of 10650 x 10~ 15 (Nishiizumi, 
2004). Measured sample “Al/^Al ratios ranged 
from 1.3 x 10 -13 to 1.7 x 10~ 12 (see GSA Data 
Repository 1 ), and AMS measurement precisions 
for 26 Al ranged from 2.4% to 12.8% (lo). We 
subtracted process blank ratios (see GSA Data 
Repository 1 ) from sample ratios and propagated 
blank uncertainties in quadrature. 


Age Calculations 

We calculated minimum limiting 10 Be and 
26 Al exposure ages (Table 2) with the CRONUS 
Earth online exposure age calculator, devel- 
opmental version 2.2, constants version 2.2 
(Balco et al., 2008), under standard atmosphere. 
We used the regionally calibrated northeastern 
North American sea-level production rates of 
3.93 ± 0.19 atoms g _1 yr 1 for 10 Be and 26.5 ± 
1.3 atoms g -1 yr 1 for 26 Al (Balco et al., 2009), 
which are thought to be the most appropriate 
production rates for west Greenland (Briner et 
al., 2012). We also used the Lal/Stone constant 
production rate model and scaling scheme, 
which defines the variation in nuclide produc- 
tion rate with latitude and atmospheric pressure 
(Lai, 1991; Stone, 2000). In CRONUS, correc- 
tions were made for latitude, elevation, sample 
thickness (ranged from 1 to 7 cm), and sample 
density (assumed 2.7 g cm -3 for granite; Table 1). 

We made no corrections for snow cover or 
till cover. Shielding by snow cover would lead 
to 10 Be or 26 Al age underestimates (Schildgen 
et al., 2005); however, seasonal snow-cover 
effects at our sample sites are likely minimal. 
Using data from a weather station in the town of 
Upernavik (www.weather-and-climate.com), 
we determined that mean temperatures are 
below freezing for 8 months of the year (Octo- 
ber through May) and that -140 mm of snow 
water equivalent falls during this time at a rate 
of -7-35 mm/month. Assuming that winter 
precipitation was added in monthly incre- 
ments, that the density of settled seasonal snow 
is 0.3 g cm -3 , and that no melting or sublima- 
tion occurred until June, shielding calcula- 
tions (Gosse and Phillips, 2001) suggest that 
reported exposure ages could underestimate 
the true age by no more than 2%. This calcula- 
tion is likely a significant overestimate because 
the areas we sampled are windswept and 
exposed during winter months and because we 
did not account for snow loss during the winter. 
However, snowfall patterns during the Pleisto- 
cene were likely different than those indicated 
by the modern weather data. 

We made no corrections for postexposure 
erosion of rock surfaces, which can lead to 
cosmogenic exposure age underestimates. 
Although boulders at low elevations did not pre- 
serve striations, the surfaces we sampled were 
glacially abraded (Fig. 2), indicating at most a 
few centimeters of erosion after the sample sites 
were exposed by deglaciation. Field observa- 
tions suggest that subaerial erosion has removed 
mass from rock surfaces at high elevations (Fig. 
3); however, because we have no reliable ero- 
sion rate estimates for these surfaces, we cannot 
correct cosmogenic ages for erosion. Thus, all 


reported ages for high- and low-elevation sur- 
faces represent minimum limits. 

Reported uncertainties reflect AMS errors 
only, which we refer to as “internal.” Because 
we are primarily interested in investigating rela- 
tionships between samples within the data set 
(i.e., a bedrock sample vs. its paired boulder), 
we use internal uncertainties in our data analy- 
sis. Additional uncertainties are introduced in 
relation to nuclide production rates, elevation 
corrections, and latitude corrections when com- 
paring the ages in this study to ages in other data 
sets. However, we can expect errors in calibra- 
tion and correction to be correlated and thus to 
affect all samples similarly because our samples 
come from a geographically limited region. 

Data Reduction and Interpretation 

Single-nuclide minimum limiting 10 Be and 
26 Al exposure ages for each sample were calcu- 
lated with the CRONUS Earth online calculator 
(Balco et al., 2008). We also considered both 
nuclides together to present minimum limiting 
exposure and burial durations (Bierman et al., 
1999). Elevation-scaled production rates for our 
samples range from 4.37 to 10.8 atoms g" 1 yrt 1 
for 10 Be and from 29.5 to 72.6 atoms g _1 yr 1 for 
26 Al. We assume a 26 Al/ 10 Be production ratio of 
6.75 (Balco et al., 2009), a 10 Be half-life of 1 .36 x 
10 6 yr (Nishiizumi et al., 2007), and a 26 A1 half- 
life of 7.05 x 1() 5 yr (Nishiizumi, 2004). Calcu- 
lations involving burial assume complete burial 
with no nuclide production. 

We compare our results to a similar study 
(Bierman et al., 1999) conducted in Baffin 
Island, Canada, because of its close proximity 
and similar latitude, elevation, and sampling 
strategy. We recalculated the production rates 
for samples presented by Bierman et al. (1999) 
in CRONUS using the northeastern North 
American production rate, and remodeled mini- 
mum limiting exposure and burial durations and 
total histories using updated decay constants so 
that the data from Upernavik and Baffin Island 
are directly comparable. 

A correction for Holocene exposure was done 
as part of the two-isotope analysis. We chose a 
value of 1 1 .3 k.y. for the Holocene exposure dura- 
tion of Upernavik samples (see “Postglacial Ice 
Retreat” section of the Discussion) and 9.5 k.y. 
for Baffin Island samples based on measure- 
ments of deglaciation age from Briner et al. 
(2009). We normalized the elevation-scaled 
production rates in CRONUS to the sea-level 
production rate, calculated the inventory of 
10 Be and 26 Al accumulated during the past 1 1.3 
and 9.5 k.y., respectively, and subtracted these 
values from the total concentrations before pro- 
ceeding with other calculations. Implementing 
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this Holocene correction shifts samples leftward 
and downward on the two-isotope plot, yielding 
shorter minimum limiting exposure durations 
and longer minimum burial durations. Although 
we used the best available data to designate a 
Holocene exposure duration, this choice intro- 
duces additional uncertainty to the two-isotope 
calculations. Using a slightly different Holocene 
exposure duration (-1-2 k.y. in either direction) 
would yield minimum limiting exposure dura- 
tions differing by several thousand years and 
minimum limiting burial durations differing by 
several tens of thousands of years. 

Retreat Rate Simulations 

To estimate the rate of latest Pleistocene ice- 
sheet retreat along the sample transect, we con- 
sidered only the youngest overlapping boulder 
sample ages ( n = 6), assuming that they inherited 
no nuclides from prior periods of exposure. We 
systematically simulated potential linear retreat 
patterns taking into account the location of each 
sample and uncertainty of each age measure- 
ment (see GSA Data Repository 1 ) in order to 
calculate a statistically most likely retreat rate. 

RESULTS 

Data Overview and Minimum Limiting 
Exposure Ages 

For the 33 samples in the data set, measured 
10 Be concentrations range from 5.0 x 10 4 to 
9.1 x 10 5 atoms g" 1 (Table 2). For the 20 sam- 
ples on which 26 Al analysis was conducted, 26 Al 
concentrations range from 4.3 x 10 5 to 5.2 x 
10 6 atoms g _1 (Table 2). Single-nuclide mini- 
mum limiting exposure ages range from 10.6 to 
104 ka for 10 Be, and from 10.9 to 89.4 ka for 26 Al 
(Fig. 7). Comparison of the two single-isotope 
concentrations for an individual sample yields 
a strong correlation (r 2 = 0.99; Fig. 8). In gen- 
eral, 10 Be and 26 Al single-nuclide concentrations 
are more similar for samples with shorter expo- 
sure times, and the 10 Be/ 26 Al ratio is close to 
the production ratio of -6.75. For longer expo- 
sure times, the concentration of 10 Be increases 
relative to the concentration of 26 Al, and the 
IO Be/ 26 Al ratio decreases below 6.75. 

The 26 Al/‘°Be ratios range from -4.3 to 7.6 
(Fig. 9). The four “control” samples used to test 
for consistency with the assumed 26 Al/ 10 Be pro- 
duction ratio of 6.75 (GU001, GU002, GU098, 
and GUI 04) fall along the constant exposure 
path (see “Using Cosmogenic Nuclides to 
Understand Glacial Exposure and Erosion” sec- 
tion) and have an average 26 Al/ 10 Be ratio of 6.9 ± 
0.2 (average ± 1 standard deviation [SD]), 
which is indistinguishable from the production 
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ratio within uncertainty. Of the remaining 16 
samples for which 26 A1 was measured, two 
(GU097 and GUI 00) have 26 Al/ 10 Be ratios 
that overlap the constant exposure path within 
uncertainty. The other 14 samples fall below the 
region of constant exposure and/or steady ero- 
sion; their 26 Al/ 10 Be ratios can be explained by 
varying lengths of burial after at least one period 
of surface exposure or by exposure under a thin 
cover of sediment or ice. 

Bedrock and Boulder Comparisons 

There is a pronounced contrast between bed- 
rock and boulder nuclide concentrations and 
calculated minimum limiting exposure ages. 
Within all 1 3 sample pairs, the bedrock sample 
has a higher concentration of 10 Be than the boul- 
der sample by more than la, considering ana- 
lytical uncertainties. Bedrock samples contain, 
on average, -2.3 times as much 10 Be as corre- 
sponding boulder samples, although this value 
varies from -1.1 to 5.1 times. Corresponding 
bedrock and boulder sample pairs have dissimi- 
lar minimum limiting 10 Be exposure ages and do 
not have a 1:1 relationship (r 2 = 0.41; Fig. 10). 
A repeated measures f-test verifies that there is 
a statistically significant difference when mini- 
mum limiting bedrock and boulder ages are 
considered in a paired comparison (p = 0.007). 
Considering bedrock and boulder samples in 


two distinct populations also results in nonover- 
lapping population distributions (Fig. 11); an 
independent samples f-test indicates that there is 
a statistically significant difference between the 
populations (p = 0.022). 

Spatial Variability 

Minimum limiting sample ages vary in the 
vertical dimension; both bedrock ( R 2 = 0.43; 
n = 20) and boulder ( R 2 = 0.4 1 ; n = 1 3) minimum 
limiting 10 Be ages increase with elevation (Fig. 
12). A one-way analysis of variance (ANOVA) 
shows that the minimum limiting bedrock 
sample 10 Be age has a weakly significant rela- 
tionship with elevation (p = 0.051), using the 
relative elevation of each sample within its own 
dipstick (Table 1). The same test performed on 
boulder samples shows that the minimum limit- 
ing 10 Be age has a more significant relationship 
with elevation (p = 0.01). Spatial variability is 
not detectable in the horizontal dimension; there 
is no relationship between minimum limiting 
sample age and distance along the 100-km-long 
sample transect. 

Five bedrock samples from the nunatak 
yielded minimum limiting 10 Be ages of 43.3- 
80.9 ka, and minimum limiting 26 A1 ages of 
35.6-67.6 ka. The lowest-elevation nunatak 
sample has the youngest minimum limiting 
single-nuclide ages (43.3 and 35.6 ka for 10 Be 


and 26 A1, respectively), while the highest- 
elevation nunatak sample has the oldest mini- 
mum limiting single-nuclide ages (80.9 and 
67.6 ka for 10 Be and 26 A1, respectively). Cal- 
culated 26 Al/ 10 Be ratios range from 5.1 to 6.2, 
overlap within uncertainty, and all fall below 
the constant exposure path. The two-isotope 
minimum limiting exposure duration is least for 
the lowest-elevation nunatak sample (43 k.y.) 
and greatest for the highest-elevation nunatak 
sample (88 k.y.). 

DISCUSSION 

Results from this study suggest that the Uper- 
navik landscape has a long and complex glacial 
history. The presence of young, Holocene-age 
boulder samples at low elevations demonstrates 
that at least some surfaces are free of cosmo- 
genic nuclides inherited from previous periods 
of exposure. However, the majority of the data 
set has older (up to ca. 100 ka) minimum limit- 
ing 10 Be and 26 A1 ages and therefore does not 
provide insight about the timing of the latest 
Pleistocene deglaciation. Using dual-isotope 
analysis and correcting for the past 11.3 k.y. 
of Holocene re-exposure, calculated minimum 
limiting exposure durations range up to 1 12 k.y., 
burial durations range up to 900 k.y., and total 
histories are as high as 990 k.y., all of which are 
indicative of a landscape that has been buried 
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Figure 7. Top: Histogram of minimum lim- 
iting 10 Be exposure ages for all 20 bedrock 
and 13 boulder samples. Bottom: Histogram 
of minimum limiting 26 A1 exposure ages for 
the 15 bedrock and five boulder samples on 
which A1 analysis was performed. 



Figure 8. Correlation of l0 Be and 26 A1 concentrations for all samples with 
dual isotope analysis ( n = 20). Error bars show la internal error; in most 
cases, error bars are smaller than the data points. The gray dashed line rep- 
resents the 1:6.75 (Be:Al) production ratio. 
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Figure 9. 26 Al/ 10 Be ratios plotted against ,0 Be concentrations for 15 bedrock and five boulder samples from Uper- 
navik, as well as seven comparable samples from Baffin Island, Canada (recalculated from Bierman et al., 1999). 
Exposure durations (in yr) are shown along the thick black “constant exposure” path. Erosion rates, in m m.y. -1 , 
are shown along the bottom of the “steady erosion” envelope. Individual erosion paths for selected erosion rates 
are shown with black dotted lines. Burial paths are shown with thin gray lines, and burial isochrones are shown 
with thin gray dashed lines. Error bars show lo internal error. 
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Figure 10. Comparison of bedrock and boulder sample ages shown by a 
correlation of minimum limiting 10 Be ages for 13 sample pairs. Error bars 
show lo internal error; in most cases, error bars are smaller than the data 
points. The gray dashed line represents a 1:1 relationship between ages. 



Figure 11. Comparison of minimum limit- 
ing bedrock and boulder sample age popu- 
lations shown by box and whisker plots 
(bedrock, n = 20; boulders n = 13). The box 
encloses the area between the first and third 
quartiles, and the horizontal line represents 
the median. Whiskers show one standard 
deviation. Samples that lie outside one stan- 
dard deviation from the mean are shown 
with an asterisk; all other samples are 
shown with dots. 


Geological Society of America Bulletin, 


Geological Society of America Bulletin, published online on 23 July 2013 as doi: 1 0. 1 1 30/B3081 3.1 

Constraining landscape history and glacial erosivity using paired cosmogenic nuclides in Upernavik, northwest Greenland 


for long periods by ice that performed little if 
any erosion. 

Glacial Period Ice Thickness and Extent 

The presence of erratic boulders on upland 
surfaces demonstrates that ice during glacial 
periods in the Upernavik region covered the 
highest peaks (-1000 m a.s.l.) at least once 
since the initiation of glaciation on Greenland, 
which occurred during the late Pliocene (Bar- 
toli et al., 2005). These glacially transported 
boulders, which do not match the local bedrock 
type (e.g., GU102, boulder, foliated granite; and 
GUI 03, bedrock, unfoliated granite), suggest 
that the Upernavik uplands were once covered 
by the Greenland Ice Sheet rather than just by 
local ice caps or ice fields. Similar observation 
of mismatched bedrock-boulder lithologies has 
been made in Melville Bugt, -300 km to the 
north (Kelly, 1980), indicating that Greenland 
Ice Sheet cover of upland surfaces was regional 
in nature. 

Deglaciation ages of the farthest outboard 
samples show that ice extended at least -100 km 
northwest of the present-day ice-sheet margin, 
providing a minimum estimate for the spatial 
extent of ice cover that was lost during the last 
deglaciation. Several of these samples came 
from low elevations (e.g., GU091 at 24 m and 
GU092 at 21 m), and their ages could repre- 
sent marine emergence rather than deglaciation. 
However, a compilation of data from around 
Greenland’s coast shows that relative sea level 
in the Upernavik region was only 0-20 m above 
modem sea level during the early Holocene fol- 
lowing deglaciation (Funder and Hansen, 1996), 


placing Upernavik in the zone of least isostatic 
rebound around Greenland. Long et al. (2006) 
found that the marine limit could be approxi- 
mated by the lowest elevation of perched boul- 
ders; since we sampled perched boulders even 
at the low elevations in Upernavik, the boulder 
samples in question are likely to have been 
above the local marine limit. Additionally, there 
is no indication in the data that the samples close 
to the marine limit are systematically younger 
than those at higher elevations. Even if these 
lowest samples were submerged, the rate of land 
surface rebound during the early Holocene was 
so rapid (Ten Brink, 1974) that the effect of ini- 
tial submergence on their exposure ages would 
be insignificant. 

Isotope Inheritance and Subglacial 
Erosion Efficiency 

This data set demonstrates that cosmogenic 
isotopes inherited from prior periods of expo- 
sure are present in most of the rock surfaces near 
Upernavik. The old minimum limiting 10 Be ages, 
positive skew of the data (Fig. 7), disagreement 
between bedrock and boulder samples (Figs. 10 
and 11), and long exposure durations modeled 
from two-isotope data (Table 2) suggest that 
most samples contain more 10 Be than could have 
accumulated during the Holocene. This excess 
10 Be is inherited from previous periods of expo- 
sure, when the climate warmed and bedrock 
surfaces and perched boulders were exposed. 
Such exposure would have occurred during 
many previous interglacial periods throughout 
the Quaternary when global ice volume was low 
(Lisiecki and Raymo, 2005). 


Upland surfaces show the longest durations of 
burial. The highest-elevation bedrock samples 
(GU096, 778 m; GU103, 998 m; and GUI 10, 
745 m) all have 26 Al/ 10 Be ratios less than -5.75, 
indicating that burial occurred at least once after 
initial exposure. The average minimum burial 
duration for samples showing burial is 324 ± 
226 k.y. (1SD, Holocene-corrected); because 
glacial periods last roughly -100 k.y., these 
long burial durations demonstrate that high- 
elevation surfaces have been repeatedly covered 
by nonerosive ice or snow. Perennial snowfields 
may explain some of this burial signal; assuming 
a settled snow density of 0.4 g cm -3 , the snow- 
pack would need to reach a depth of -10 m to 
effectively shield rock surfaces from spallation- 
induced nuclide production. Complete burial 
by cold-based glacial ice or thick snowfields, 
or partial burial by thin snowfields or till, could 
have occurred following initial exposure (Fig. 
5). However, the presence of discordant bedrock 
and boulder ages indicates that the boulders 
were emplaced at some point after the initial 
exposure and burial of the bedrock; thus, the 
most likely explanation for mismatched mini- 
mum limiting bedrock-boulder ages is burial by 
cold-based ice at least once after initial exposure 
(Fig. 5). 

The presence of inherited 10 Be indicates that 
subglacial erosion rates were low, and that gla- 
cial ice was not capable of removing all 10 Be 
inherited from preceding interglacial periods 
(Fig. 5). Furthermore, the relationship between 
minimum limiting sample age and elevation 
(Fig. 12) suggests that glacial ice removed 
more material from low elevations than it did 
from high elevations (Briner et al., 2006). This 
interpretation is consistent with field observa- 
tions. Low-elevation areas display streamlined 
landforms indicative of glacial erosion (Fig. 2), 
while high-elevation areas preserve evidence of 
extended surface weathering (Fig. 3). Elevation- 
related erosion intensity is likely caused by the 
relationship between ice thickness and subgla- 
cial temperatures; low-elevation land surfaces 
were covered by thicker ice, at or near the pres- 
sure melting point, which would have been 
more erosive than the thinner cold-based ice that 
covered high-elevation surfaces. 

While elevation-related erosion intensity 
governs the overall trend of sample ages with 
elevation, there is still noise in the data. The 
sample with the longest modeled total his- 
tory (GUI 10, 989 ka) is from an elevation of 
745 m a.s.l.; there are numerous higher-elevation 
samples with shorter total histories. The five 
samples from the nunatak (dipstick 7; Fig. 6) 
do not show a correlation between elevation and 
modeled total history, maybe because all of the 
samples have similar elevations (within 100 m), 



Elevation (m a.s.l.) 

Figure 12. Correlation of minimum limiting sample age and elevation for 
bedrock samples (n = 20) and boulder samples ( n = 13). Error bars show lo 
internal error; in most cases, error bars are smaller than the data points. 
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and ice erosivity did not vary appreciably over 
this thickness. Deviations from the trend of 
increasing sample age with elevation likely arise 
from local variability in glacial erosive power, 
differences of subaerial erosion during inter- 
glacial periods, rock hardness, small shielding 
effects, and the possibility that ice lingered lon- 
ger at lower elevations. 

Some boulder samples also preserve inherited 
10 Be, although not as much as bedrock samples. 
Boulders were quarried farther inland and then 
transported coastward to their present locations. 
They may have originated from inheritance- 
bearing outcrops, which were originally covered 
by nonerosive ice but subsequently covered by 
erosive ice capable of quarrying. They also may 
have been plucked from outcrops far inland 
from the coast, where ice is thicker and more 
likely warm-based, before being transported 
into a zone of cold-based ice around the mar- 
gin. There is appreciable spread in the minimum 
limiting 10 Be ages of boulder samples (ranging 
fourfold, from 10.6 to 45.7 ka), possibly docu- 
menting a variety of different sources, erosion 
mechanisms, or transport histories. 

The low efficiency of subglacial erosion we 
infer in Upemavik is consistent with findings in 
other high-latitude areas. Inheritance of cosmo- 
genic nuclides from prior periods of exposure 
has been documented in Greenland (Goehring 
et al., 2010; Hakansson et al., 2008; Kelly et 
al., 2008), the Canadian Arctic (Bierman et 
al., 1999; Briner et al., 2003; Marquette et al., 
2004), Scandinavia (Goehring et al., 2008; Har- 
bor et al., 2006; Stroeven et al., 2002b), and Ant- 
arctica (Lilly et al., 2010; Sugden et al., 2005). 
These regions are characterized by cold mean 
annual temperatures that would favor the exis- 
tence of an ice sheet partially frozen to the bed. 

However, not all high-latitude regions have 
cold-based ice. The presence of inherited 
nuclides in Upemavik differs notably from the 
Jakobshavn fsbrae region, 500 km south. There, 
no inheritance of cosmogenic isotopes from 
prior periods of exposure was detected, regard- 
less of sample type or elevation (Corbett et al., 
201 1). The lack of inheritance there may reflect 
thick ice in Jakobshavn Isbrte, which eroded 
sufficient material from bedrock outcrops to 
remove any inherited isotopes. The pronounced 
contrast between inferred subglacial erosion 
efficiency in Upemavik (low) and Jakobshavn 
(high) indicates that landscapes buried by ice 
in or around large ice streams likely experience 
more efficient erosion than those buried by thin- 
ner, more slowly flowing ice. Climate probably 
also plays a role, since the mean annual tem- 
perature in Upemavik is ~2.5 °C colder than 
in llulissat (www.weather-and-climate.com). 
Subglacial erosion efficiency, therefore, may be 


controlled by local factors such as ice thickness 
and ice-flow velocity, as well as larger-scale 
regional climatic conditions. 

Landscape Evolution 

Parts of the Upemavik landscape have been 
preserved through multiple glacial-interglacial 
cycles. High-elevation bedrock surfaces have 
minimum limiting total histories of hundreds 
of thousands of years, suggesting that the high- 
lands were minimally altered through repeated 
burial by glacial ice. Burial dominates the land- 
scape history; most high-elevation surfaces have 
minimum limiting burial durations of -300- 
500 k.y. and minimum limiting total histories 
of -500-600 k.y., suggesting that these surfaces 
experienced burial with limited erosion -75% of 
the time. 

The high-elevation surfaces on Baffin Island, 
-600 km to the west, have longer minimum 
limiting exposure durations, burial durations, 
and total histories than similar surfaces in 
Upemavik (using recalculated data from Bier- 
man et al., 1999; Figs. 9 and 13). On average, 
bedrock surfaces on Baffin Island record total 
histories almost two times as long (-900 k.y.) 
as those in Upemavik. However, despite dif- 
ferences in duration, the relative role of burial 
appears to be similar between Baffin Island and 
Upemavik. The Baffin Island samples have a 
burial to total history ratio of 0.79 ± 0.04 (aver- 
age ± 1 SD); the Upemavik samples (taking 
into account only bedrock samples at elevations 
>600 m in order to match the Baffin Island data 
set) have a burial to total history ratio of 0.83 ± 
0.07 (average ± 1 SD). Thus, despite recording 
different durations of Quaternary history, both 
upland surfaces have experienced burial for the 
same portion of time. The value of -0.8 is likely 
a reflection of the fraction of time represented 
by glacial periods and expanded ice volume 
throughout the latter part of the Quaternary. 

In Upemavik, the presence of cosmogenic 
nuclide inheritance in bedrock samples deep 
within fjords raises the question of how the 
fjords formed, if not by glacial erosion. Given 
that glaciers can form a recognizable U-shaped 
valley from a preexisting V-shaped valley in 
as little as 10,000 yr (Harbor, 1992), even one 
or two glacial periods of thick, highly erosive 
ice could cause significant reshaping of the 
landscape. However, it is probable that ice 
cover in Upemavik has been cold-based for 
the majority of the Quaternary, and has not 
been a major player in landscape develop- 
ment. More likely, the origin of these fjords 
is tectonic, and thus glacial activity has been 
only a recent overprint upon an older structure 
(England, 1987; Glasser and Ghiglione, 2009; 


Swiff et al., 2008). Bonow et al. (2007) sug- 
gested that uplift of the west Greenland coast 
occurred during the late Miocene and Pliocene, 
as seafloor spreading was initiated in the North 
Atlantic. Continued erosion created planar sur- 
faces, which were uplifted to create the low- 
relief highlands visible in western Greenland 
today. Fluvial erosion incised these surfaces as 
rivers responded to changing base levels, creat- 
ing a landscape of planar highlands dissected 
by deep valleys (Bonow et al., 2007). Repeated 
burial by cold-based, nonerosive glacial ice 
during the Quaternary has likely done little 
to modify this landscape; this implies that the 
fjords in Upemavik may be relict features, pre- 
served from the Miocene or Pliocene (Bonow 
et al., 2007). 
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Figure 13. Histograms of minimum limiting 
exposure durations, burial durations, and 
total histories from highlands in Upemavik, 
Greenland (this study), and Baffin Island, 
Canada (recalculated from Bierman et al., 
1999). The Upemavik data include only 
bedrock samples from elevations >600 m in 
order to be comparable to the samples from 
Baffin Island. All samples have been cor- 
rected for 11.3 k.y. (Upemavik) or 9.5 k.y. 
(Baffin) of Holocene exposure and scaled 
for elevation-based differences in produc- 
tion rates. 
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Postglacial Ice Retreat 

Determining the timing of ice retreat in the 
deep fjords after the last glacial period is pos- 
sible because a small subset of the samples 
appear to have simple exposure histories. The 
six youngest boulder samples have overlapping 
10 Be ages of 1 1.3 ± 0.5 ka, and we interpret this 
to be the best estimate for the deglaciation of 
the exposed land surface near Upernavik. This 
reported age is the mean exposure age ± one 
standard deviation; systematic uncertainties 
in production rate would shift the central ten- 
dency slightly but not change the relationship 
between samples. This cosmogenic estimate for 
deglaciation is consistent with other estimates of 
minimum limiting deglaciation age from radio- 
carbon dating discussed in the “Study Site and 
Previous Work” section. 

The low-elevation boulder samples show no 
trend in age along the 100-km-long sample tran- 
sect, suggesting that significant lateral retreat of 
the ice margin (-100 km) likely occurred within 
several hundred years. The modeled statistical 
maximum likelihood ice retreat rate was -120 m 
yr 1 , with possible rates as low as 80 m yr 1 and 
as high as 240 m yr 1 (1 SD). Rapid retreat of 
the ice margin was likely facilitated by calving, 
since ice in the Upernavik region would have 
been largely marine-based due to local topog- 
raphy, which consists of small islands separated 
by large fjords (Fig. 6). These results support the 
proposition that areas of marine-based ice can 
undergo periods of rapid collapse in a warming 
climate (England, 1976; Mangerud et al., 2013; 
Straneo et al., 2012). 

The retreat of the ice-sheet margin in Uper- 
navik shows a pattern of retreat similar to that of 
the Helheim Glacier in southeastern Greenland. 
There, Hughes et al. (2012) used 10 Be exposure 
dating on bedrock and boulder samples and 
detennined that samples collected along the 
entire length of the -60 km fjord had indistin- 
guishable exposure ages within two SD. The 
mean age of boulders indicates that retreat took 
place very quickly at 10.8 ± 0.3 ka. The timing 
of rapid ice loss from Helheim Fjord (Hughes 
et al.. 2012) overlaps with the estimate of rapid 
ice loss from Upernavik presented here (1 1.3 ± 
0.5 ka) and suggests that the episode of abrupt 
early Holocene retreat seen in Upernavik is not 
spatially restricted. The rapid episode of glacial 
retreat documented in both Upernavik and in 
Helheim Fjord appears to have occurred earlier 
than retreat in Ilulissat, 500 km to the south, 
where the present-day coastline did not become 
ice free until ca. 10.3 k.y. B.R (Corbett et al., 
2011; Young et al., 2011). 

The rate of ice-margin retreat in Upernavik is 
similar to some but higher than many estimates of 


retreat rates for large bodies of ice from the same 
time period (see Fig. 9 in Corbett et al., 2011). 
In Ilulissat, 500 km to the south, calculated rates 
of ice-margin retreat were -100-110 m yr 1 
through Jakobshavn Isfjord, its tributary fjord, 
and its adjacent bay (Corbett et al., 2011; Long 
and Roberts, 2003; Long et al., 2006; Young et 
al., 2011). In Sisimiut Fjord, central western 
Greenland, retreat rates were -20-50 m yr 1 
(Rinterknecht et al., 2009). In North America, 
retreat rates of -60 m yr 1 were calculated for 
Sam Ford Fjord on Baffin Island (Briner et al., 
2009). Ice retreat rate, therefore, appears to be 
controlled by local and regional factors. 

Rapid ice loss from the Upernavik area ca. 
1 1 .3 ± 0.5 ka may have been driven by external 
climate forcing (Fig. 14). Oxygen isotope ratios 
from the Greenland Ice Sheet Project 2 (GISP2) 
ice core (Stuiver and Grootes, 2000) and from 
the Greenland Ice Core Project (GRIP) and 
Dye3 ices cores (Dahl-Jensen et al., 1998) 
show that this time period was characterized by 
rapid warming. This warming may have coin- 
cided with the end of the Younger Dryas cold 
event (Alley, 2000; Taylor et al., 1993), which 
is thought to have occurred rapidly in the Arctic 
over several abrupt subdecadal pulses at 1 1 ,645 ± 
200 and 11,612 ± 200 yr B.P. (Taylor et al., 
1997). However, any causal relationship is spec- 
ulative because only six of the samples appear 
to record deglaciation age, and because the aver- 


age deglaciation age of the Upernavik fjords has 
an uncertainty of 500 yr (1 SD). 

CONCLUSIONS 

Analysis of 10 Be and 26 Al in paired bedrock 
and boulder samples near Upernavik, northwest 
Greenland, reveals a complex landscape history. 
Bedrock and boulder samples have discordant 
10 Be and 26 Al minimum limiting ages, suggest- 
ing that nuclides inherited from previous peri- 
ods of exposure are present. High-elevation sur- 
faces have old single-nuclide minimum limiting 
exposure ages and record at least several hun- 
dred thousand years of burial after initial expo- 
sure. These surfaces, which exhibit pronounced 
surface weathering features, have likely been 
preserved under nonerosive, cold-based ice or 
perennial snowfields over the course of numer- 
ous glacial cycles. The existence of these relict 
surfaces suggests that, despite repeated gla- 
ciation, certain high-latitude landscapes have 
remained virtually unchanged for hundreds of 
thousands or even millions of years. Conversely, 
low-elevation surfaces have younger minimum 
limiting ages, although bedrock samples still 
preserve inheritance and evidence of burial. 
Boulders at low elevations indicate that degla- 
ciation along the 100-km-long sample transect 
took place rapidly at a rate of -120 m yr 1 at ca. 
1 1.3 ka, supporting the idea that large bodies of 
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Figure 14. Deglaciation timing and paleoclimate context. Thin black lines 
show the probability density functions for each individual boulder sample 
(the oldest boulder sample, 44.3 ka, has been omitted from this plot for bet- 
ter visibility). The thick black line shows the summed cumulative probabil- 
ity for the entire boulder data set. The dark gray line represents § ls O values 
from the Greenland Ice Sheet Project 2 (GISP2) ice core from Greenland’s 
summit; these data, available online at http://depts.washington.edu/qil/ 
datasets/, are courtesy of M. Stuiver and P.M. Grootes at University of 
Washington. Minimum limiting "'Be ages are plotted against the left axis, 
while GISP2 5 ls O values are plotted against the right axis. 
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ice can disintegrate on time scales that would be 
noticeable and meaningful to humans. 
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